Non-Gaussian fluctuations near the QCD critical point 
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We study the effect of the QCD criticaf point on non-Gaussian moments (cumulants) of fluctu- 
ations of experimental observables in heavy-ion collisions. We find that these moments are very 
sensitive to the proximity of the critical point, as measured by the magnitude of the correlation 
length £. For example, the cubic central moment of multiplicity {(5N) 3 } ~ £ 4 ' 5 and the quartic cu- 
mulant {(6N) 4 ) C ~ £ 7 . We estimate the magnitude of critical point contributions to non-Gaussian 
fluctuations of pion and proton multiplicities. 
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Introduction — Mapping the QCD phase diagram is 
one of the fundamental goals of heavy-ion collision ex- 
periments. QCD critical point is a distinct singular fea- 
ture of the phase diagram, the existence of which is a 
ubiquitous property of QCD models []J based on chiral 
dynamics. Locating the point using first-principle lattice 
calculations is a formidable challenge and, while recent 
progress and results are encouraging, much work needs 
to be done to understand and constrain systematic er- 
rors 0. If the critical point is located in the region ac- 
cessible to heavy-ion collision experiments it can be dis- 
covered experimentally. The search for the critical point 
is planned at the Relativistic Heavy Ion Collider (RHIC) 
at BNL, the Super Proton Synchrotron (SPS) at CERN 
and the future Facility for Antiproton and Ion Research 
(FAIR) at GSI. 

This Letter focuses on the experimental observables 
needed to locate the critical point in heavy- ion collisions. 
Locating the point requires a scan of the phase diagram, 
by varying the initial collision energy y/s. The charac- 
teristic signature is the non-monotonous behavior, as a 
function of \fs, of the experimental observables sensitive 
to the proximity of the critical point to the point where 
freeze-out occurs for a given y/s. d, 0] 

The most characteristic feature of a critical point is 
increase and divergence of fluctuations. Most fluctuation 
measures discussed to-date can be related to quadratic 
variances of event-by-event observables, such as particle 
multiplicities, net charge, baryon number, particle ratios, 
or mean transverse momentum in the event [5]. Typi- 
cally, the singular contribution to these variances induced 
by the proximity of the critical point is proportional to 
approximately £ 2 , where £ is the correlation length which, 
in the idealized thermodynamic limit, would diverge at 
the critical point [J] . The magnitude of £ is limited triv- 
ially by the system size, but most stringently by the 
finite-time effects due to critical slowing down [4j, |6| . The 
observation Q that the correlation length may reach at 
most the value of 2 — 3 fm, compared to its "natural" 
value of 1 fm, may make discovering the critical non- 
monotonous contribution to such fluctuation measures a 
challenging task, if the measures depend on £ too weakly. 

In this Letter we point out that higher, non-Gaussian, 
moments of the fluctuations are significantly more sensi- 



tive to the proximity of the critical point than the com- 
monly employed measures based on quadratic moments. 
We explore this observation quantitatively and calculate 
the critical point contribution to selected fluctuation ob- 
servables. 

Illustration — The point can be illustrated using a 
description of the fluctuations based on the probability 
distribution of an order parameter field - any field which, 
by quantum numbers, can mix with the critical mode a 
- the mode developing infinite correlation length at the 
critical point. 

We choose the maximum of the probability distribu- 
tion to be at a — 0. The probability distribution P[a] 
can be written as 



P[a] ~exp{-fi[a]/T} : 



(1) 



where f2 is the effective action (free energy) functional 
for the field a, which can be expanded in powers of a as 
well as in the gradients: 



O = 



2 V ' 2 3 4 



(2) 

Near the critical point m a <C T , so the mode a can be 
treated as a classical field. 

Calculating 2-point correlator (a(x)a(O)) we find that 
the correlation length £ — m^ 1 . For correlation functions 
of the zero momentum mode ctq = / d 3 x o~{x)/V we find 
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= y[2(A 3 2 -A 4 k 8 



The critical point is characterized by £ — > 00. The central 
observation in this Letter is that the higher moments 
(cumulants) K3 and K4 diverge with £ much faster than 
the quadratic moment «2- uM 

Of course, the fluctuations of the critical mode are 
not measured directly in heavy-ion collision experiments. 
These fluctuations do, however, influence fluctuations of 
multiplicities, momentum distributions, ratios, etc. of ob- 
served particles, such as pions or protons, to which the 
critical mode couples 0] • The purpose of this Letter is 
to determine the magnitude of these effects. 
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FIG. 1: Diagrammatic representation of the contribution to 
the three-particle correlator from the critical mode a. Wavy 
lines represent propagators of the a field, each contributing 
factor l/m&, crossed circles represent insertions of Sn p into 
the correlator Eq. ((4]) - see Ref.[3] for details. 



FIG. 2: Diagrammatic representation of the critical mode 
contribution to the connected four-particle correlator. Same 
notations as in Fig. [1] 



given by 



Critical contribution to experimental observables - 
We shall now estimate the effect of the critical point 
fluctuations on the observables such as the pion multi- 
plicity fluctuations. Using similar approach, it should be 
straightforward to construct corresponding estimates for 
such observables as charge, proton number, transverse 
momentum fluctuations, etc., as well as to take into ac- 
count acceptance cuts. 

We shall focus on the most singular contribution, pro- 
portional to a power of the correlation length £. This 
contribution can be found using an intuitive picture de- 
scribed in Ref. [4j. In this picture one considers a joint 
probability distribution for the occupation numbers n p 
of observed particles (such as pions) together with the 
value of the the critical mode field a (more precisely, its 
zero-momentum mode (To), the latter treated as classical. 
Due to coupling of the critical mode of the type amr, the 
fluctuations of the occupation numbers receive additional 
contribution, proportional to the corresponding correla- 
tion functions (moments) of the fluctuations of o~q given 
by Eq. (|3|). In this Letter, however, it will be more conve- 
nient to use instead more formal diagrammatic method 
developed in Ref. [7(. 

Cubic cumulant — The 3-particle correlator receives 
the following most singular contribution from the a fluc- 
tuations, given by the diagram in Fig. [TJ 



(Sn pi 5n P2 5n T 



2A 3 
V 2 T 



(4) 



11)7, 



UJ pi UJ. 



P2 W P3 



where, subscript a indicates that only contribution of 
the critical mode is considered and, as in Refs. d 0> 
we denoted airit coupling by G and introduced a short- 
hand notation for the variance of the occupation number 
distribution: v p = n p (l ± n p ), where the "+" is for the 
Bose particles. 

Since the total multiplicity is just the sum of all occu- 
pation numbers and thus 



5N = ^2 5n p> 



(5) 



the cubic moment of the pion multiplicity distribution is 



((6N) S 



pi J p2 J p3 



(6) 



where J = J d 3 p/(27r) 3 . Since ((SN) 3 } scales as V 1 it is 
convenient to normalize it by the mean total multiplicity 
N which scales similarly. Thus we define 



N 



and find 



2A 3 G 3 



T ml 



(7) 



(8) 



Quartic cumulant — The leading contribution to the 
connected 4-particle correlator is given by the sum of 
two types of diagrams in Fig. [2j 



(Sn pl Sn P2 5n P3 Sn 
G 
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G\ 4 v 2 pi v 2 p2 v 2 p3 v 2 pi 



(9) 



Pi ^Vl ^P3 W P4 



The quartic cumulant of multiplicity fluctuations is 
given by 



{(6N) 4 } C = V 4 (5n pi Sn P2 5n P3 6n P4 ) c (10) 

J P1P2P3P4 

This cumulant also scales as V 1 in thermodynamic limit. 
As in Eq. ([7]) we define a ratio whose V — > oo limit is 
finite: uj 4 {N) = ((SN) 4 ) C /N, and find 



X 2 

2 7T — X4 




(11) 

Estimate of the effect — In order to estimate the mag- 
nitude of the effect, we need to estimate the values of 
the coupling constants A3 and A4. The main uncertainty 
in the estimate will come, however, from the uncertainty 
of the value of G, which enters in a large power. This 
constant is known only roughly - the estimate was made 
in Ref. [J]. Therefore, a crude estimate for A's suffices. 
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FIG. 3: Illustration of the possible relative position of the 
critical point and the locations of the freeze-out points for 
different values of the initial collision energy \/~s in a heavy- 
ion collision. 



Near the critical point both A3 and A4 vanish with a 
power of £ given by (neglecting 77 <C 1): 

\ 3 = \ 3 T-(T0~ 3/2 , and A 4 = A 4 ■ (TO -1 , (12) 

where dimensionless couplings A3 and A4 are universal, 
and for the Ising universality class they have been mea- 
sured (see, e.g., Ref.|8] for a review). A3 varies from 
to about 8 depending on the direction of approach to 
the critical point (crossover or first-order transition side) . 
The coupling A4 varies from about 4 to about 20. Since 
the freeze-out occurs somewhere between these two ex- 
tremes (as illustrated in Fig. [3|) we shall pick some mid- 
range values for our estimates. The main point is the 
strong dependence of the effect on £. 

Putting together estimates of A, G (from [4|) and £ 
(from 3) we find for pions at T w 120 MeV (in full 
acceptance, for a single pion species, see also Eq. (|20"|) ): 



w 4 (iV w ) c 



G 



300 MeV 



3 fm 



9/2 



(13) 




Because of large powers of the coupling G, which is known 
only poorly, the uncertainty in this result is significantly 
larger than that of similar estimates of the quadratic mo- 
ments of fluctuations [1] . 

The most significant feature of this result is the strong 
dependence on £ which makes the cubic and quartic cu- 
mulants very sensitive signatures of the critical point. 

Another example: proton multiplicity fluctuations - 
The above analysis carries over to the proton multiplicity 
fluctuations. The fluctuations of the net proton number 
is a good proxy to the baryon number fluctuations, whose 



magnitude, proportional to the baryon number suscepti- 
bility, must diverge at the critical point Ref. . This sus- 
ceptibility, as well as "kurtosis", K4/K2, have been stud- 
ied on the lattice using Taylor expansion around [ig = 
and in QCD models Here, for simplicity, we shall 

present the results for the proton only multiplicity, which 
is easier to measure in experiments. 

To adapt equation such as Eq. Q to protons, one needs 
to substitute G with the coupling g of the critical mode 
a to protons (times the mass m p of the proton): G — ► 
gm p . The estimate for this coupling can be taken from 
the sigma-model to be roughly g w m p / w 10. The 
variance of the occupation number distribution is — 
n p (l - n p ) and fi p = ( e ( w p~^ B )/ T + l) -1 , where \ib is 
baryochemical potential. There is also a factor of 2 k ~ 1 for 
Wfe (as in Eq. ([20]) ) because of the proton spin degeneracy. 
Putting this together one finds, e.g., for protons at SPS 
freeze-out conditions (T,fi B ) w (168,266) MeV [Ti| 



w 3 (iV P ) 



Ui{N p ) 




(15) 



(16) 



Note that the effect is much larger on the proton multi- 
plicity fluctuations, compared to thepion multiplicity. 

Similar to quadratic fluctuations [9j], the exponents in 
Eqs. ([T5jl . (fTtj)) agree (up to 77 <§C 1) with the critical 
exponents of the baryon number cumulants dictated by 
scaling and universality: 



((SN £ 



VT 



-,-x d k P(T, hb) 
8(1% 



fc(5-?7)/2-3 



(17) 



Mean transverse momentum — From the expression 
for the correlators Eq. (j4|) or (J9j) one can similarly esti- 
mate the effect of the critical point on other observables, 
for example, higher moments of the fluctuation of mean 
transverse momentum px- For example, the cubic mo- 
ment k 3 {8pt) of the mean px distribution around the 
all-event mean pr can be expressed as 



k 3 {Spt) = ({pi 



Pt) 



E 

Pl,p2,P3 



([Pi]t ~Pt) 



(18) 



x Qpzlr ~ Pt){[p 3 ]t ~ PT)(8n pi Sn P2 Sn P3 ) 



and estimated using Eq. ((4]). Normalizing this variable 
similarly to the variable F proposed in Ref. [3| removes 
N scaling and makes it less sensitive to the effect of the 
flow ("blue shift" of momenta): 



k 3 (pt) 
NvUc(pt) 



(19) 



where Vi nc (pT) is the variance of the inclusive (single- 
particle) pt distribution. We leave this to future work. 
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Comments and discussion — It is worth noting that, 
even though the £ dependence of L04 is stronger, its 
measurement involves subtraction of two contributions, 
((6N) 4 )— 3((5N) 2 ) 2 , each of which is order N times larger 
than their difference, which might dilute the signal-to- 
noise ratio in experimental measurement. 

Since the freeze-out occurs, generically, somewhat off 
the crossover line, as illustrated in Fig. 02 one should 
expect the critical point contribution to fluctuations to 
be skewed. In this case, the deviations from the Gaussian 
shape are dominated by the cubic moment, or W3. 

What is the sign of W3(iV) CT ? One can anticipate it by 
using the following, admittedly crude, argument. The 
skewness of the distribution of the order parameter near 
the critical point is a "shadow" of a second peak. This 
peak corresponds to the phase on the other side of the 
first-order transition line (quark-gluon plasma phase at 
higher T and [is on Fig. [3]). This phase has higher 
entropy and baryon number, thus fluctuations of these 
quantities must be skewed toward higher values: 103 > 0. 
Since pion and proton numbers are rough proxies to en- 
tropy and baryon number respectively their skewness 
should be also positive. 

What is the "natural" , background value one should 
expect for, e.g., cj 3 ? For a gas of classical free parti- 
cles (Poisson distribution) ujs(N) = 1. Bose statistics 
increases this by a factor ^(A^be = (1 + n p )(l + 2n p ), 
e.g., approximately 1.3 for pions at T = 120 MeV. 

More importantly, quantum statistics only corre- 
lates fluctuations of particles of the same species, thus 
U3(N)be is the same for all charge (i.e., N = N n + +N^- ) 
and single charge multiplicity fluctuations. In contrast, 
the critical point contribution correlates also ir + with 7r~, 
thus making W3 4 times larger for all charge vs single 
charge (for A\+ = N v - ) : 



u k (N v+ +N v -) a = 2 k - 1 u> k {N„+) 



(20) 



Eq. ([201) can help separate critical point contribution 
from contributions due to quantum statistics. 

It is important to note that other sources may and do 
contribute to the skewness and kurtosis: remnants of ini- 
tial fluctuations, flow, jets - to name just a few obvious 
contributors. Identification and evaluation of these con- 
tributions is a task well beyond this paper. This serves 
to emphasize that the energy scan of the QCD phase 
diagram is needed to separate such background contri- 
butions from the genuine critical point effect, the lat- 
ter being non-monotonous function of the initial colli- 
sion energy y/s as the critical point is approached and 
then passed. The fact that non-Gaussian moments have 
stronger dependence on £ than, e.g., quadratic moments, 
makes those higher moments more sensitive signatures of 
the critical point. 
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